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ABSTRACT 

Aims. We present the detection, identification and calibration of extended sources in the deepest X-ray dataset to date, the extended 
Chandra Deep Field South (ECDF-S). 

Methods. Ultra-deep observations of ECDF-S with Chandra and XMM-Newton enable a search for extended X-ray emission down 
to an unprecedented flux of 2 x 10“*® ergs s“* cm“^. By using simulations and comparing them with the Chandra and XMM data, we 
show that it is feasible to probe extended sources of this flux level, which is 10,000 times fainter than the first X-ray group catalogs of 
the ROSAT all sky survey. Extensive spectroscopic surveys at the VLT and Magellan have been completed, providing spectroscopic 
identification of galaxy groups to high redshifts. Furthermore, available HST imaging enables a weak-lensing calibration of the group 
masses. 

Results. We present the search for the extended emission on spatial scales of 32" in both Chandra and XMM data, covering 0.3 square 
degrees and model the extended emission on scales of arcminutes. We present a catalog of 46 spectroscopically identified groups, 
reaching a redshift of 1.6. We show that the statistical properties of ECDF-S, such as logN-logS and X-ray luminosity function are 
broadly consistent with LCDM, with the exception that dn/dz/dD test reveals that a redshift range of 0.2 < z < 0.5 in ECDF-S is 
sparsely populated. The lack of nearby structure, however, makes studies of high-redshift groups particularly easier both in X-rays 
and lensing, due to a lower level of clustered foreground. We present one and two point statistics of the galaxy groups as well as 
weak-lensing analysis to show that the detected low-luminosity systems are indeed low-mass systems. We verify the applicability 
of the scaling relations between the X-ray luminosity and the total mass of the group, derived for the COSMOS survey to lower 
masses and higher redshifts probed by ECDF-S by means of stacked weak lensing and clustering analysis, constraining any possible 
departures to be within 30% in mass. 

Conclusions. Ultra-deep X-ray surveys uniquely probe the low-mass galaxy groups across a broad range of redshifts. These groups 
constitute the most common environment for galaxy evolution. Together with the exquisite data set available in the best studied part 
of the Universe, the ECDF-S group catalog presented here has an exceptional legacy value. 

Key words, galaxy groups - galaxy evolution 
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1. Introduction 


Detection of extended X-ray emission is an important source of 
information on the hot intergalactic medium of groups and clus¬ 
ters of galaxies. A sample of X-ray groups recovered by deep 
surveys is a unique resource to improve our understanding of 
low-mass groups as well as distant clusters. It also provides in¬ 
formation on the common environment of massive galaxies. 


The advent of Chandra and XMM-Newton has ele¬ 
vated galaxy group research to a new level, with large 
catalogs of X-ray selected grou p s now availa b le for many 
surveys jFinoguenov et al.l 2007L 2010L_ 2009I: George,e^^ 


2(^ : lAdami et al.ll20 1 It IConnellv et al.ll20 1 21: lErfanianfar et al 


20131) . The first studies using those catalogs have already re¬ 
vealed substantial differences in the galaxy population of galaxy 
groups: compared to galaxy clusters, gro ups haye more baryons 
locked in galaxies ([Giodini et al.| |2 009 | ). and haye more star¬ 
forming galaxies (iGiodini et al. l2012l : IPoDesso et al.ll2012h . The 
redshift eyolution of the star-formation rate in groups has been 
found to differ fro m clusters, approach ing the field leyel at inter¬ 
mediate redshifts (iPonesso et al.ll201^ . Diyersity of the optical 
propert ies of high-z groups has been reported by iTanaka et al.l 
(I2013bl) . 


The ability of X-rays to characterise galaxy groups in terms 
of their mass and yirial radius enables a robust separa t ion of 
mass and radial trends in galaxy formation. IZinaro et al.l (l2014t) 
showed that a fundamental difference exists between X-ray de¬ 
tected groups and group-like density regions, where enyiron- 
mental processes related to a massiye dark matter halo are more 
efficient in quenching galaxy star formation with respect to 
purely density related processes. In particular, the rapid eyo¬ 
lution of galaxies in groups with respect to group-like density 
regions and the field highlights the leading role of X-ray de¬ 
tected groups in the cosmic quenching of star formation. Use 
of groups proyides a direct estimate of the halo occupation dis¬ 
tribution, which are not affected by the sample yariance, as well 


(Smolcic et al.l2011 

; George et al.ll201 ll 

2012UFeauthaud et al.l 

1201 2l Allevato et al. 

120121 Oh etalJl2014 

). 
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Fig.l : Combined Chandra and XMM exposure map of ECDF-S area. 
The contours represent levels of 0.1, 1, 2, 4, 8 and 12 Ms effec¬ 
tive Chandra ACIS-I on-axis exposure. 


X-ray galaxy groups, however, have proven to be more diffi¬ 
cult objects to study at X-rays, compared to clusters. Therefore, 
the role of surveys in finding galaxy groups is particularly 
unique. The large depths required to study the galaxy groups are 
rewarded by their high volume abundance. One has literally just 
to stare at any direction for sufficiently long time to find them. 

Among all X-ray surveys, the Extended Chandra Deep Field 
South (ECDF-S) is by far the deepest X-ray survey on the 
sky. The galaxy group catalog recovered in this work is there- 
fore of unique impor tance. Following the pioneering work of 
iGiacconi et al.l (l2002l) . this paper presents a systematic account¬ 
ing of the extended X-ray emission in the ECDF-S area, based 
on a factor of 10 deeper data, with an equivalent Chandra ACIS- 
I exposure of 16 Ms in the central (CDF-S) area (see ^for de¬ 
tails). 

This paper is structured as follows: in ^we describe the X- 
ray analysis; in Owe describe the identification of X-ray galaxy 
groups; in Owe present the modelling of the X-ray detection of 
galaxy groups; in Owe discuss the properties of the groups and 
present the one-point statistics; in Owe present the clustering 
analysis and our modelling of the bias; in ^4.11 we present the 
modelling of the observed emission in the entire ECDF-S field, 
based on the identification of groups and their properties; in O 
we present the stacked weak lensing profile; in Owe discuss the 
ECDF-S superstructure at a redshift of 1.6. Results are discussed 
in section 


2. Data and analysis technique 

2.1. XMM-Newton and Chandra data reduction 


The ECDF-S area has been a frequent target of X-ray ob¬ 
servations with both Chandra and XMM . After the first IMs 
Chandra observation (IGiacconi et al.ll2()()^ . the area was named 
the Chandra De ep Field South. The extension of the CDF-S 
surve y to 2 Ms ( Fuo et ani2008l) and later to 4Ms of exposure 
time (fXue et al.ll201 ll) . via a large Director’s Discretionary Time 
project, has now provided our most sensitive 0.5-8 keV view of 
the distant AGNs and galaxies. This paper does not include the 
3Ms Chandra observations of the field taken in 2014. 

For the detection of extended sources, a dominant con¬ 
tribution to t he sensitivity is p rovided by ultra-deep XMM 
observations (iRanalli et al.ll20ll) . obtained under several pro¬ 
grams, most importantly a 3Ms Very Farge Program (PI: Andrea 
Comastri). For the XMM data analysis we have followed 
the prescription outlined in iFinoguenov et al.l (l2007l) on data 
sc reening and backgr ound evaluation, with updates described 
in iBielbv et al.l (l2010l) . After cleaning those observations from 
flares, the resulting net total observing time with X MM-Newton 
are 1. 946Ms for the pn (for a description see IStriider et al.l 
(|20^ ). 2.55 2MS for MOSl, and 2.530Ms for MOS2 (for a de¬ 
scription see iTurner et al.l (120011) '). For detecting the extended 
emission on arcminute scales, the sensitivity of each MOS is 
similar to Chandra ACIS-I, while pn detector is 3.6 times more 
sensitive. We adopt the Chandra ACIS-I units of exposure, 
adding XMM EPIC pn exposures with a weight factor of 3.6. 
We refer to it as an effective Chandra exposure, as it corresponds 


* All observed values quoted through this paper, are calculated 
adopting a A CDM cosmological model, with //„ = 70 km s“* Mpc^', 
Q.M = 0.24, Q.\ = 0.76 (but see the modelling for testing Planck cosmo¬ 
logical parameters. We quote all X-ray fluxes in the [0.5-2] keV hand 
and rest-frame luminosities in the [0.1-2.4] keV band and provide the 
conhdence intervals on the 68% level. FK5 coordinates used through¬ 
out. 
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Fig. 2: Sensitivity of Chandra and XMM towards the detection of X-ray emission on 32" scales. Contours show the levels of 1.2,2,3, and 5x10“'* 
ergs s“' cm“^ and provide the intensity scale in the image. Left: all Chandra observations. The sensitivity does not reach the deepest contour. 
Middle: All XMM observations. Right: Chandra plus XMM. 


to the time required by Chandra to achieve the same sensitivity 
on > 32" scales. In FiglUwe show the resulting exposure map 
of the survey. The peak exposure of the survey is 16 Ms. 

In the Chandra analysis we apply a conservative event 
screening and modelling of the quiescent background. We fil¬ 
ter the event light-curve using the lc_clean tool in order to 
remove normally undetected particle flares. The background 
model maps have been ev aluated with the prescription of 
iHickox & Markevitchl (l2006h . We estimated the particle back¬ 
ground by using the ACIS stowed position observation^ and 
rescaling them by the ratio . The cosmic background 

flux has been evaluated, by subtracting the particle background 
maps from the real data and masking the area occupied by 
the detected sources. The rapid changes in the Chandra point 
spread function (PSF) as a function of off-axis angle produces a 
large gradient in the resolved fraction of the cosmic background, 
which is the primary source of systematics in our background 
subtraction. 

For cataloguing the groups, we also include the ECDF-S 
data (iLehmer et al.ll2005h . which consist of four Chandra ACIS- 
I pointings, 250ksec each, defining the square shape of the ex¬ 
posure and sensitivity maps in Figs [1] and |2] However, a simple 
addition of the Chandra ECDF-S and Chandra CDF-S data re¬ 
sults in a reduction in the quality of background subtraction in 
the CDF-S area, coming from the outer part of ECDF-S ACIS-I 
data. So for the final analysis we include the dataset with re¬ 
moved ECDF-S ACIS-I data in overlap with the CDF-S ACIS-I 
data and use the simulations of the held (§4), which reproduce 
the low sensitivity of the corners of ACIS-I ECDF-S mosaic. 

2.2. Point source subtraction 

To detect and study faint extended sources, we must begin with 
the removal of flux pro duced by the point sources, following 
iFinoguenov et al.l(l2009l) . We model the position-dependent PSF 
of each instrument and subtract the model from the XMM, the 
Chandra CDF-S and the Chandra ECDF-S mosaics separately, 
using the flux map of point sources derived from each mosaic. 
In subtracting the point sources, we operate with a flux distri¬ 
bution on small scales, as reconstructed using wavelets, without 
an attempt to catalog the sources or to use existing point source 

^ http://cxc.cfa.harvard.edu/contrib/maxim/acisbg 


catalogs. The point source emission is resolved in Chandra, but 
can be confused for XMM. For XMM we remove the flux from 
point sources down to flux levels of 10“'® ergs s“* cm“^ in the 
0.5-2 keV band, below the a corresponding confusion limit for 
XMM (10“'® ergs s“' cm-^). 

For Chandra, the point source contribution to the spatial 
scales in excess of 16 arcseconds in the central (within 3' radius 
from the average aim point) detector area is negligible, while the 
ratio of flux on scales of 8-16 arcsecond to over 16 arcsecond 
can be approximated as a co nstant at large (> 3') off-axis an- 
gles IFinoguenov et alJ (l2009l) . The point source subtraction pro¬ 
cedure separates out the flux below 8 arcseconds and uses the re¬ 
maining flux detected within the 8-16 arcsecond scale to predict 
the residual contamination on scales above 16 arcseconds. The 
systematic effects associated with variation in the flux attributed 
to a given scale by wavelets (noisy sources have less flux de¬ 
tected on smaller sc ales) were mitig a ted by using the calibrated 
wavelet program of IVikhlinin et al.l (Il998h and applying three 
levels of flux reconstruction, with 4, 30 and 100 sigma detection 
thresholds and using different flux scaling for each significance 
level. We have verified that our flux maps for Chandra contain 
a contribu t ion fro m all ~ 750 catalogued AGNs and galaxies in 
IXue et all (1201 Ih . The residuals due to asymmetric PSF shapes 
were quantified and added as a systematical error to preclude 
their detection. For XMM, the selectio n of spatial scales used for 
point source flux has been explained in iFinoguenov et all (1201 Oh 
and consists in absence of off-axis behaviour in the encompassed 
flux ratios below and above 16". The effect of subtracting off 
the contribution from point sources is extremely important for 
XMM, as illustrated in Fig[3 The number of extended features 
is reduced by a factor of 2, and the appearance of an XMM im¬ 
age on 32 arcsecond scales becomes similar to that of Chandra. 

In Figl^we show the sensitivity towards the detection of ex¬ 
tended emission after the contamination from both background 
and point sources have been removed. In the 0.5-2 keV band, the 
Chandra data alone reach fluxes of 2 x 10“*® ergs s“* cm“^, while 
XMM data alone reach 1.2 x 10“'® ergs s“* cm“^. The combined 
dataset reaches similar depths to that of the XMM alone, but over 
larger area. The quoted flux corresponds to the detection cell of 
0.7 arcmin^. Detailed simulations of the detection are discussed 
in §4. 

We performed an analysis on simulated XMM maps of point 
sources, presented in iBrunner et alJ (l2008l) for similarly large 
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Fig. 3: Wavelet reconstruction of the XMM image on 32" - 128" scales without (left panel) and with (right panel) the flux removal coming from 
the wings of point sources. The number of apparent sources changes by a factor of 2. 


XMM exposures in the Lockman Hole, detecting no extended 
emission in the simulated maps containing the detected by XMM 
point sources. The higher sensitivity of Chandra towards the 
detection of point sources allows us to make a statistical as¬ 
sessment of the effect of sub-threshold (for XMM) AGNs to¬ 
ward the detection of extended emission. Performance of XMM 
observations was accompanied by deepening the Chandra data 
within one year from each other, which makes Chandra maps 
suitable for XMM point source contamina tion analysis, limit¬ 
ing the effect of AG N long-term variability (ISalvato et al.ll20lil 
iPaolillo et alJl200^ . We have computed the variation of unre¬ 
solved point source flux on the detection scales for XMM, us¬ 
ing the Chandra image, masking out the sources detected in the 
XMM analysis. The constructed Chandra flux map has been fur¬ 
ther smoothed with a Gaussian of 16" width, approximating the 
effects of the XMM PSF. In the map, the uniform distribution 
of the faintest point sources results in nearly constant emis¬ 
sion, which we subtract following the procedure for local cos¬ 
mic background estimates for XMM, while bright sources and 
clustered sources make an enhancement. We find the contami¬ 
nation by point sources unresolved by XMM to the flux of iden¬ 
tified extended sources is below the 5% level of the extended 
source’s flux. The highest peaks in the contamination map are 
associated with stand-alone sources near the (XMM) detection 
threshold, which by chance happened not to coincide with any 
of the detected groups and would contribute 30% to the faintest 
group flux. The importance of these so urces is even higher in 
shallow surveys (iMirkazemi et al.ir2015h . to a degree requiring 
matched detection thresholds between point-like and extended 
sources, effectively removing faint extended sources from con¬ 
sideration. The importance of point source remo val in XMM 
data is mentioned also i n other cluster publications (iHilton et alJ 
lIoTolilPierre et^lIoTl . 

Our procedure for point source removal has been exten¬ 
sively tested on the real observations and is tuned for the actual 
XMM PSF. We have previously tested our pipeline on the simu¬ 


lation s of the Lockman Hole (iHenrv et al.ll2010l : iBrunner et al.l 
l2008h . finding no residuals. For the ECDF-S, we can extend 
those tests to an image a factor of 5 deeper and include the effects 
of sub-threshold AGNs down to fluxes of 10“*^ ergs s“* cm“^, 
based on the deep Chandra catalogs. In FiglHwe show the simu¬ 
lated image and the residuals detected on 32" - 128" scales. We 
have simulated point sources flux and the background for each of 
the XMM pointings, and followed the procedure for background 
and point source subtraction. 

A total of 16 extended sources have been detected in the 
simulated 0.3 square degree mosaic image, while only point 
sources were used as an input. These fake extended sources 
correspond to large-scale distribution of unresolved sources by 
XMM and each source is made of a combination of typically 
7 AGNs inside the source and lack of AGNs on either part of 
the source. We also performed a detection of simulated point 
sources, adding an error associated with the extra flux due to the 
extended sources. The number of detected fake sources has not 
decreased substantially (15), 8 of those are in the CDF-S area. 

Finally, since the positions of the simulated sources are real, 
they should correspond to an actual extended source in XMM. 
The number of such detections in XMM mosaic is 3. This is due 
to the fact that most of the fake sources being close to the flux 
limit of 2 X 10“'® ergs s“* cm“^, where detection is affected by 
the confusion on extended emission. The 3 detected fake sources 
have a flux of 2, 3, 5 xl0“*® ergs s“' cm“^, with a correspond¬ 
ing flux error of 1.2 x 10“'® ergs s“' cm“^, which agrees with 
the detected flux by XMM at those positions. None of these fake 
sources were identified as galaxy groups and entered the final 
catalog. However, they have contributed to a reduction in the 
identification rate by 6%. In Fig|5]we overlay the contours of 
detected extended emission over the simulated point source con¬ 
tamination image. 

In Figl6]we show the signal-to-noise ratio obtained for the 
final joint dataset (excluding ECDF-S Chandra data) after sub¬ 
tracting the background and detected point sources. The white 
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Fig. 4: Simulated XMM mosaic image of point source and background 
emission in the ECDF-S. 




Fig. 5: Simulated residuals after the point source removal and back¬ 
ground subtraction. Contours show the detected X-ray emission. 
We identify three detected false sources (highlighted by dashed 
circles). None of these sources were identified as a galaxy group. 


part of the image corresponds to zero or negative signal. The 
grey and black parts of the image correspond to an area with 
significant flux, which occupies a substantial (20%) part of the 
image. There are three large sources; one in the east, associ¬ 
ated with a nearby group; one in the north, associated with a 
nearby cluster having a peak outside the area of ECDF-S, but 
seen clearly in the ACIS-S chip that was on during the obser¬ 
vation; the third source, which is near the center, is due to con- 



-1 0.6 7.f i.t &.« A.6 10 I? li 15 


Fig. 6: Signal-to-noise of the XMM data after point source removal and 
smoothing with a 16" Gaussian kernel. The color bar shows the 
correspondence between the color and the significance of the 
emission, starting with white for -Icr. 


fusion of several groups with overlapping virial radii. We will 
return to the modelling of the image in ^4.11 

2.3. Source extraction 

The sensitivity of the source detection depends critically on the 
background per resolution element. The level of the background 
per unit area is comparable between Chandra and XMM. On 
small scales, the XMM PSF leads to large corrections for the en¬ 
compassed flux of the source, which reduces the effective XMM 
sensitivity towards point sources. On scales selected for the anal¬ 
ysis in this paper, the PSF does not affect the source flux, but 
there is an induced background due to a distribution of AGN 
counts by larger PSF of XMM. These differences support a con¬ 
sideration of separate Chandra and XMM searches for the ex¬ 
tended sources, in addition to a joint search. 

Sources found in deep X-ray surveys ar e primarily AGNs 
and distant galaxies (iBrandt & Hasinged2005h . Groups and clus- 
ters of galaxies only ac count for 10% of the cosmic flux (e.g. 
iFinoguenov et al.ll2007h . Their emission on arcminute scales re¬ 
quires different detection methods versus compact sources. Most 
techniques to date refer to detection of galaxy groups and clus¬ 
ters as extended sources. 

The term extended emission is however loosely defined. To 
some extent any astrophysical emission results from objects that 
are not singular and so it is only a question of how extended 
the emission is. Emission on scales of a few arcseconds in the 
survey data appears to stem from the cores of the groups. X-ray 
jets, galaxy mergers and even individual galaxies. 

An important characteristic of group X-ray emission is a 
correlation between its intensity and angular extent; the emission 
typically covers a sizable fraction of the /? 5 oo radius that can be 
derived based on the observed flux and a known source redshift. 
Groups of galaxies that are sufficiently bright to be detected in 
X-rays, exhibit emission on arcminute scales even at the high¬ 
est redshifts accessible to the deepest surveys like the ECDF-S. 
As the detection is background limited, and given the shape of 
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Fig. 7: Left: XMM detection of extended emission on a 32 arcsecond scale. Right: Chandra detection of extended emission on a 32 arcsecond 
scale. Contours, which are the same in both panels, show the extended emission detected in the combined Chandra and XMM images on 
the 32 and 64 arcsecond scales. Full ECDF-S field of 0.3 sq.degrees is shown. 


the surface brightness profile of galaxy groups, the emission on 
smaller scales is more easily detectable. The adoption of spatial 
scales of 32" is therefore a trade-off between signal-to-noise on 
one hand and both telescope characteristics and source identi¬ 
fication, on the other. The depths of the ECDF-S preclude us¬ 
ing large spatial scales, as due to the high number of extended 
sources the emission is confused on the arcminute scales. 

In Fig 1?] we compare the final detection map with the indi¬ 
vidual maps obtained by Chandra and XMM. The most signif¬ 
icant sources appear in both maps. For the final detection, we 
combined the residual maps of the Chandra ECDF-S, CDF-S 
and XMM ECDF-S. The practical issue of the combining maps 
with different pixel sizes is handled using the TERAPIX S WARP 
software. We co-add the residual counts without any weight, co¬ 
add the exposure maps re-normalised to differences in the effec¬ 
tive areas of the instruments and add the error maps in quadra¬ 
ture. The sensitivities of Chandra and XMM towards the X-ray 
emission in the 0.5-2 keV band also depend on the spectrum of 
the group emission, while in adding the data we can only assume 
a typical ratio of the sensitivities. Large differences in the ratio 
of sensitivities occur only if the emission is primarily at energies 
below 0.7 keV, where also the differences between pn and MOS 
are large. In Fig0we compare the XMM and Chandra fluxes 
for the sample. We use the effective exposure units, in which the 
count-rate of XMM and Chandra are similar. We view Figl8]as a 
characterisation of the scatter introduced by our attempt to merge 
XMM and Chandra raw counts, which is of the order of 0.2 dex. 
A few bright objects are located at the outskirts of the observa¬ 
tions and also occupy a large area, leading to instrument-specific 
differences in the background prediction. 

3. Identification of galaxy groups 

All sources in our catalog are X-ray selected, using the emission 
from outskirts of the groups, typically exceeding 100 kpc scales 
(with any exception from this criteria duly noted), uniquely iden¬ 
tifying galaxy groups even at low luminosities. X-ray data alone 
are not sufficient for source identification and thus our effec- 



Fig. 8: Flux comparison between Chandra and XMM within the area 
covered by the 4Ms Chandra CDF-S. The solid line shows the 
1:1 correspondence. 38 extended sources with significant flux 
measurement in both Chandra and XMM data are shown. The 
errors on XMM fluxes are similar to the plotted Chandra errors 
and are omitted from the plot for clarity. 


tive survey sensitivity is a combin ation of both X-ray and op- 
tical/NIR sensitivities. For examnle. lBielbv et al.l (1201 Oh demon¬ 
strated that deep NIR data are essential to identify distant groups 
and clusters of galaxies. In this work, we combine the ultra-deep 
X-ray observations of the ECDF-S with the exquisite optical- 
nearlR photometric and spectroscopic data available in the field . 

We run our r ed sequence finder (iFinoguenov et al.l 120101 : 
iBielbv et al.ll2010t) around the central portions of all the X-ray 
group candidates. We base our red sequence sear ch on the Penn 
State photometric redshift catalog described in iRaffertv et al.l 
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(1201 ih . We first extract galaxies around a redshift of interest 
by applying \Zphot - z\ < 0.1. We then cou nt galaxies around 
the mo del red seq uence constructed w ith the iBruzual & CharloB 
(l2003h model ('see lLidman et al.l|2008l for details). When count¬ 
ing, we use a Gaussian weight in the form of 



/ colon^ohs - color„,odeiiz) 

\ ^i,obs 


magi,„bs - 

^ mag 

expj-(^)2j, (1) 

where coloriobs and mag jobs are the color and the magnitude of 
the i-th observed galaxy, a-j^ots is the observed color error in 
colori^obs, color„iodei{z) is the model red sequence color at the 
magnitude of the observed galaxy, mag*^^^^^j(z) is the character¬ 
istic magnitude based on the model, which is tuned to reproduce 
the observed characteristic magnitudes, cr,nag is the smoothing 
parameter and is set to 2.0 mag, r, is the distance from the X- 
ray center and cr,. is another smoothing para meter with 0.5 Mpc. 
In our earlier work (iFinoguenov et alJ2010h . we adopted cr^ = 1 
Mpc, but here we apply a smaller window of 0.5 Mpc because we 
search for both smaller and more abundant (we therefore need to 
reduce the chance association) systems. The significance of the 
red sequence around an X-ray source is computed with respect 
to the mean and variance of the number of red galaxies measured 
at random positions in the same field. 

Since different colors are sensitive to red galaxies at dif¬ 
ferent redshifts, we adopt the combination of colors and mag¬ 
nitudes summarised in Table [T] We use t he publicly availabl e 
MUSYC photometry in the ECDF-S area (iGawiser et al.ll200^ . 
which is slightly smaller than the full X-ray coverage. In the 
GOOD S area, we use the deeper public catalog from the MUSIC 
survey dGrazian et alJl20d^lSantini et alJl2009l) . Our experience 
shows that we need to go down to ~ M* -F 1 to securely identify 
a red sequence. High redshift systems lack faint red galaxies, 
but the red sequenc e is often seen down to that magnitude (e.g., 
iTanaka et alJlioO'^ . The MUSYC data for ECDFS is not deep 
enough to identify z > 1.5, and thus high-z identifications are 
not yet complete at present. The MUSIC data is deep enough to 
see systems at z = 2 and beyond. In fact, we have identified two 
z ~ 1.6 groups as discussed below. 

One may worry that a red sequence finder introduces a bias; 
it may miss groups dominated by blue galaxies. But, we note that 
a red sequence finder misses only groups in which the red frac¬ 
tion is significantly smaller than the field. Suppose the red frac¬ 
tion in a group is the same as the field, a group is an over-density 
of galaxies by definition and thus there is a larger number of red 
galaxies within a small volume, which will then be detected by 
a red sequence finder. It is an interesting question if groups with 
a lower red fraction than the field exist at high redshifts. They 
may, but recent observations of z ~ 2 systems, es pecially those 
in the p rocess of forming, show red sequence ('e.g.. lTanaka et alJ 
l2013bl) . This might indicate that red sequence is a ubiquitous 
feature of groups and clusters since an early epoch. 

In the identification process, we have made an extensive 
use of spectroscopic redshifts available in the field. The X- 
ray data used in this work covers a 0.3 square degree area, 
which is larger than the 0.1 square degree area of GOODS-S, 
where extensive public ESO spectroscopic surveys have been 
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Fig. 9: Distribution of core radii of the detected simulated ex¬ 
tended sources. The shades show the detected sources in the 
"Confusion” run, with shades of grey illustrating the overlap 
of sources. The dashed and solid lines show the location of 
90% and 50% detection completeness level. The core radii are 
uniquely determi ned by the mass and reds hift of the halo, using 
the tabulations of iFinoguenov et aP ( 120071) . 


carried out (e.g. iBalestra et al.ll201(ih . Two spectroscopic cam¬ 
paigns have been used to remedy this situation: ECDF-S follow¬ 
up th rough devoted ESO and Keck efforts (ISilverman et al.l 
1201 oh . and since 2009, the fol low-up of groups h as been car¬ 
ried out by the ACES project (ICooper et al.ll201^ . which is a 
large program on the Magellan telescope. We have complied 
a spectroscopic catalog with a high sampling rate (60% down 
to i=22 AB mag) from these efforts. We replace red sequence 
redshifts with spectroscopic redshifts where available. Large 
amount of spectroscopy available in the field, enables a search 
for the spectroscopic galaxy groups, with most massive ones 
having a good correspondence to th e location of X-ray emission 
(iDehghan & Johnston-Holfit3l2014h . 

4. Modelling of the X-ray source detection 
procedure 

In this section we provide the validation of the X-ray detection 
method. Readers not interested in the technical details of the X- 
ray detections may skip to section ^ 

Our method of detection of extended sources differs from 
other X-ray surveys. Most X-ray cluster surveys aim to fit a 
symmetrical beta model with a fixed beta to a list of extended 
source candidates, resulting in the determination of the cluster 
core radius. This modelling of the surface brightness profile is 
later used to infer the total flux of the cluster within some radii 
(iPacaud et alj|200'^ iLlovd-Davies et~aDl201 ll) . If we revisit the 
origin of the method, the reasoning for using it comes primar¬ 
ily from low-redshifts, where the core is well resolved, while 
outskirts of the cluster are not observed. The high reliance of 
X-ray surveys on the core properties of clusters has been ar¬ 
gued by a number of studies as a weakness, as it introduces a 
large scatter in cluster selection, favouring the detection of clus¬ 
ters with strong cool cores . On the other hand, it has been ar¬ 
gued (IVikhlinin et al.ll2()()^ that cluster outskirts exhibit a much 
smaller scatter with total mass, as w itnessed in the low-scatter 
of core-excised Lx (lMaughanll2007h . The reported low-scatter 
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Fig. 10: Distribution of the beta parameter of the detected simulated ex¬ 
tended sources. The shades show the detected sources in the 
’’Confusion” run, with shades of grey illustrating the overlap 
of sources. The dashed and solid lines show the location of 
90% and 50% detection completeness level. The value of beta 
is uniquely determined by the mass and redshi ft of the halo, 
using the tabulations of iFinoguenov et ^ l l2007b . 


measurements were obtained using a simple aperture flux. It 
therefore seems logical to pursue a method of cluster detection, 
which would only be sensitive to flux coming from outskirts. 
In addition, resolving the core of a high-z group is only possi¬ 
ble with the on-axis PSF of Chandra, while for XMM resolving 
cores below 10" is both incomplete (iLlovd-Davies et al]l201 ih 
and is subject to contamination from point sources (iPacaud et al l 
l2007h . For redshifts above 0.5, detecting the outskirts beyond 
half of the R^qq value is typical. 

The low scatter of core-excised Lx suggests that the clus¬ 
ter surface brightness profile can be modelled as a sum of two 
profiles, one describing the core and the other describing the 
outskirts with the ratio of core to Rsqq radii and values of beta 
similar to that of merging clusters. While this is yet to be ver¬ 
ified, it implies a low-scatter scaling between the detected flux 
in the fixed aperture and the core-excised Lx in the region en¬ 
closing the flux calculation. Deviations from this assumption 
would violate the published low scatter of gas mas s presented 
for a wide range of overdensities, from 2500 to 500 dAllen et alJ 
l2008t IVikhlinin et al.ll200^ lOkabe et al.ll2010t) . which brackets 
the overdensities important for this work. Given these consider¬ 
ations, we assume that we can restore the core-excised flux of 
the group based on the detection of group outskirts. 

In order to model the X-ray detection of groups in ECDF- 
S, we explore different aspects of the group detection. Using 
the sensitivity map, presented in Fig|2l we calculate the limit¬ 
ing group mass that can be detected for each area of equal sen¬ 
sitivity, using a grid of redshifts and adopted scaling relations 
with total mass. For each limiting mass, redshift, and volume, 
corresponding to equal sensitivity areas and steps of the redshift 
grid, we generate simulated groups with masses according to the 
mass function, defined by LCDM with the Plan ck cosmolog¬ 
ical parameters (iPlanck Collaboration et al.|l2014ll . In adopting 
the set of parameters we select the Planck CMB only constraints 
(no BAG), which are also close to WMAP9. The values of the 
cosmological parameters assumed are Q,„ = 1 - Qa = 0.30, 
erg = 0.81,/j = 0.7. 


With about 100 groups expected from the cosmology, we 
would not be able to sample well all the parameter space im¬ 
portant for detection and so the expected numbers were boosted 
by a factor of 100, constrained by the time required to perform 
the simulations. The positions of the groups were randomised 
within the sensitivity area and the redshifts were randomised 
within the resolution of the redshift grid. In simulating the halos 
with masses in excess of 10'"^M©, we do not follow the shape of 
the mass function in detail, but calculate the integral of the mass 
function above 10and upon boosting and randomising the 
total number of simulated systems, and we assign the 10 
mass to all such sources. This creates an upper boundary in the 
point distribution visible in Figs.9-13. 

As a second step, we run direct simulations of the group de¬ 
tection. For each gro up, we used the total mas s to establish R^qq 
and the tabulations of iFinoguenov et akl (l2007h to predict the pa¬ 
rameters of the beta model. The limiting flux of the detection 
is translated to the limiting value of beta and core radii. Figs|9] 
and [To] illustrate the values of core radii and beta as a function 
of redshift for detected sources. We show the curves of incom¬ 
pleteness, calculated for the CDFS area, showing where we start 
loosing sources, which is determined by the flux on the detec¬ 
tion scales. Most high-redshift groups have an expected extent 
of their X-ray emission (/? 5 oo) comparable to the detection scale 
used, while the values of their core radii can only be resolved 
by Chandra. As a result of the point source removal applied, the 
cores of the simulated groups are removed as well and the detec¬ 
tion is only sensitive to the flux at the group outskirts, while the 
extent of the simulated detection approaches Rsqq- At z < 0.3 the 
core of the group becomes detectable, and variations in the in¬ 
ner group surface brightness become important. The simulations 
generate the group profile, and projects it on the exposure map. 
This provides a model to further pixel-wise randomisation of a 
number of photons detected and the model for errors which we 
add to the survey noise map. 

Tab|2 summarises the results of source detection simula¬ 
tions for two choices of scaling relations, the COSMOS one and 
adopting a 30% higher mass for given total L^, to mimic the 
effect of our calibration uncertainties. The later are termed as 
’’Scaling” runs. The basic run is termed ’’Sensitivity only”. We 
consider the additional effects of confusion and confusion-t-PSF. 
We cull our input group list for simulation near detection bound¬ 
ary, based on the experience with sampling the parameter space. 

To simulate the effect of confusion, we reran the simula¬ 
tions, adding simulated groups to the actual ECDF-S image. We 
require the peak of the emission to be within 16" of the original 
center. We remove the area within 16" from the peaks on the 
X-ray image, as those would always be detected. These peaks 
occupy 5% of the total area. As seen in Tab|2| (“confusion” run), 
accounting for confusion leads to further reduction in the total 
number of sources, however due to the statistical nature of the 
detections and enhanced detection of the emission at places near 
the existing sources produces a small number (2%) of detections 
in this run, which were not obtained in the previous one. 

To simulate the differences in the source detection between 
Chandra and XMM, we performed another round of simulations 
with confusion, in which we convolved the source profile with 
the XMM PSE and in the detection procedure we introduce a 
step of flux removal from large scales, based on the detection 
on small scales. We find (Tab|2| ’’Confusion-tPSE” run) that the 
main result is a 2% increase in the detection rate. Thus, we con¬ 
clude that XMM PSE only marginally inhibits the detection of 
groups in our algorithm, compared to running it on the Chandra 
data. The origin of the effect is due to tuning of the flux removal 
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Fig. 11: Flux-redshift plane of the ECDF-S sample (filled circles with 
error bars). The grey shades show the distribution of the param¬ 
eters of the detected groups in simulations, with Planc k cosmo¬ 
logical parameters i Planck_Collal^ationet al .11 20141) and the 
scaling relations of Leauthaud et alJ <2010h used. Due to the 
limited spatial scales used, at low redshifts the effective sensi¬ 
tivity towards the total flux is lower. The upper boundary on the 
flux distribution shows the combination of the ECDF-S survey 
volume and cosmology. The dashed and solid lines show the 
location of 90% and 50% detection completeness level. 



z 


Fig. 12: Luminosity-redshift sampling in the ECDF-S. Grey shadow¬ 
ing indicate the density of detected groups in simulations, with 
Planc k cosmological parameters l Planck_Colh^oration et alJ 
l2014ll and the scaling relations of lLeauth'auTet^J ( 120101) used. 
Solid circles show the parameters of ECDF-S groups. The 
dashed and solid lines show the location of 90% and 50% de¬ 
tection completeness level. 


for the point sources, leaving in a fraction of the flux from the 
group cores scattered by XMM PSF to outskirts of the groups. 
In the ’’Scaling” runs, we only considered the effects of change 
in the scaling and confusion (Tabl2]i. 

We can compare the results of the simulations also to the 
combined XMM-i-Chandra catalog. We use the results of the 
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z 

Fig. 13: Mass-redshift sampling in the ECDF-S. Definition of mass is 
done with respect to the mean density and is scaled by the 
Hubble constant. Grey shadowing indicate the density of de¬ 
tected grou£sJn_simulutionvw^ Planc k cosmological param¬ 
eters l lPlanckCollaboration^e^ allI20l4l and the scaling rela¬ 
tions of lLeauthaiMet^J l l2010h used. Solid circles show the pa¬ 
rameters of ECDF-S groups. The dashed and solid lines show 
the location of 90% and 50% detection completeness level. 


’’Confusion” run. We illustrate the detections in FigHT] The sim¬ 
ulations allow us to show the expected completeness of X-ray 
group detection as a function of luminosity or a group mass, 
which we illustrate in FigsfT2land[T3] We show the complete¬ 
ness curves calculated for the CDFS area in Figsl9lfT3] In Tab[3 
we summarise for the three representative flux levels the proper¬ 
ties of the survey in terms of contamination and present the es¬ 
timates of the completeness at z=0.6, with account for the effect 
of confusion. In FiglTT]we also see a gradual loss of sensitivity 
towards the group detection with increasing redshift at z > 0.6 
caused by the reduction in the angular size of R 50 Q. 

Outside the radius of three core radii, for a given slope of 
the surface brightness profile, the scaling of the emission from 
one spatial bin to another does not depend on the actual value 
of the core radius, as can be shown analytically. Since we con¬ 
sider the variation of the central luminosity of the X-ray group 
as a source of scatter, ignoring this variation shall be understood 
as the low-scatter part of Lx, just-like the core-excised Lx- The 
absolute value of the Lx can deviate even from the average Lx 
for groups of a given flux and redshift. For the purpose of infer¬ 
ring the group mass, this requires calibration, for which we use 
external methods, such as clustering and weak lensing. Thus, 
even if the actual group parameters would systematically devi¬ 
ate from the assumed ones and exhibit the scatter, we can still 
rely on our method of assigning the total mass. The actual scal¬ 
ing relation will however be method (and thus instrument) de¬ 
pendent. In our modelling, we go from the cosmology to the 
mass function, and then to the expected Lx given the calibra¬ 
tions suitable for our parametrisation of the surface brightness 
profile and later evaluate the detection. Should we change the 
surface brightness profile parametrisation, the scaling relations 
would have to be changed to compensate for the change in the 
flux in the detection cell. The actual variation of the flux on large 
scales for a given mass is expected to be as small as the reported 
behaviour of core-excised Lx, which is 7% for clusters, accord- 


9 























A. Finoguenov et al.: Ultra-deep catalog of X-ray groups in the Extended Chandra Deep Field South 



Fig. 14: Wavelet reconstruction of the simulated image of extended 
emission in the ECDF-S area on scales of 0.5-2 arcminutes. 
The contours of the observed X-ray emission on matched spa¬ 
tial scales are overlaid in black. Contours not aligned with sim¬ 
ulated X-ray emission correspond to unidentified sources. The 
contours do not show the largest scales of the emission, but 
similarity between the model image and the signal-to-noise im¬ 
age in Fig.6 is clear. 

ing to iMaughar] (l2007h . which is negligible, compared to the 
statistical scatter for the simulated (and used) 4cr detection limit. 

Our modelling is performed under an assumption of no 
evolution of the fraction of surface brightness associated with 
0.1f?5oo. Existing statements in the literature, indicate that if any, 
the cool core contribution to the total flux is reduced. Thus we 
believe that our assumptions are conservative. For comparison 
with literature, we note that the importance of the emission in¬ 
side the core radius is much higher for steep beta values, like 0.6, 
which is typically assumed for and is a characteristic of massive 
clusters. 

4.1. Understanding the effect of group outskirts in expiaining 
the X-ray image 

The high spatial density of sources, identified in the ECDF- 
S exposures, should result in largely overlapping emission on 
large scales. To test this effect, we use the identified systems to 
model the X-ray image on large spatial scales. We assumed a 
beta model for each of the groups with core radii equal to 10% 
of the virial radius and a slope p - 0.6. High (3 - 0.6 values as¬ 
sumed, can be viewed as conservative for estimating the source 
confusion on large scales, as the surface bright profile for each 
source drops fast. The normalisation is chosen to match the aper¬ 
ture flux of the source. The simulated exposure approximately 
matches the achieved sensitivity. A flat exposure map and 5" 
PSF are adopted, using the SIXTE (ISchmid et alJl2010h Athena 
WEI set-up. These differences are not important for making our 
point. To compare the simulated image with the observed one, 
we applied the same wavelet reconstruction procedure and in 
FigOcompare the detected emission on 0.5-2 arcminute scales. 
The revealed similarity in the image is quite striking. The details 
of the arcminute-scale variation in the X-ray emission are well 


reproduced. This emission caused problems for estimating the 
sky background, in the northern and eastern part of the survey, 
leading us to use the central vs western part of the survey for in¬ 
field estimates of instrumental vs sky background components. 
The complex bright structures on 2' scales in the ECDF-S are re¬ 
produced as an effect of confusion on large scales. And even the 
complex appearance of the sources on arcminute scales seems to 
be sufficiently modelled as the confusion of several sources (e.g. 
the ’’Fudge” source is a combination of four galaxy groups). 

5. Galaxy groups in the ECDF-S field 

5.1. The Group Catalog 

In this section we describe our catalog of 46 X-ray galaxy groups 
detected in the ECDF-S field as well as estimates for the 5 com¬ 
ponents of the Kurk structure (iKurk et al.l l2009h . In the cata¬ 
log (TablUi we provide the source identification number (col¬ 
umn 1), lAU name (column 2), R.A. and Deck of the X-ray 
source in Equinox J2000.0 (3^), and redshift (5). The clus¬ 
ter flux in the 0.5-2 keV band is listed in column (6) with the 
corresponding 1 sigma errors. The flux has units of 10 ergs 
cm~ ^ s~* and is extrapolate d to an iteratively determined Rsqq 
(see lFinoguenov et al.ll2007l for details). The aperture determin¬ 
ing the flux has been defined by the shape of the emission on 
32" scales, unless it has been manually redefined to avoid con¬ 
tamination from other extended sources (cases where this is not 
possible have flag=4). The total net XMM-i-Chandra counts in 
the flux extraction region are given in (7). The rest-frame lumi¬ 
nosity in the 0.1-2.4 keV band in units of 10"^^ ergs s“^ is given 
in (8), where the K-correction assumes the temp e rature from 
the scaling relations adopted in iFinoguenov et al.l (l2007h . The 
choice of the energy ba nd is driven by the avai lable calibrations 
of the Lx - M relation (iFeauthaud et al.ll2010l) . yielding (col.9) 
an estimated total mass, M 200 , defined with respect to the critical 
density, with only the statistical errors quoted. Sy stematic errors 
due to scatter in the scaling relations are ~ 20% (lAllevato et al.l 
l2012h and the uncertainty on the calibration is 30%, as discussed 
in ^ and ^ The corresponding R 200 in arcminutes is given in 
column (10). Column (11) lists the source flag and the number of 
spectroscopic member galaxies inside R 200 , used to evaluate the 
mean spectroscopic redshift, is given in column (12). In Column 
(13) w e provide the predicte d galaxy velocity dispersion based 
on the lCarlberg et alJ d 19971) virial relation using our total mass 
estimates. A comparison be tween these and actual m easured val¬ 
ues of Vdi'ip is presented in lErfanianfar et aP (l2014l) . The errors 
provided on the derived properties are only statistical and do not 
include the intrinsic scatter in the Lx-M relation and the system- 
atics associated with the extrapolation of the scaling relations to 
lower luminosities at similar redshifts. In ^we successfully ver¬ 
ify these masses by means of a clustering analysis to a precision 
below the 0.2dex uncertainty of individual mass estimates due to 
the scatter in M - Lx relation. In ^we also successfully verify 
the mass calibration by means of stacked weak lensing analysis. 

While a number of groups we report on were previously dis¬ 
covered by Chandra, their emission has only been probed out to 
much smaller radii, and so it was much more uncertain as a char¬ 
acterisation of the group properties. This poses a trade-off for 
optimising future telescope performance, as detection benefits 
from high angular resolution, while the characterisation benefits 
from low background and collecting area. 

There is an issue related to the definition of the extended 
source flux, correspo ndin g to quotation of th e source flux. In 
iGiacconi et al.l (120021) and iBauer et al.l (l2004h the detected flux 
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is quoted, while in iFinoguenov et al.l (l2007i 1201 Oli the full flux 
of the source is quoted. These can be different by a factor of a 
few. Using the large spatial scales, one reduces the amount of 
extrapolation on the flux and therefore removes a large separa¬ 
tion between the ob served and referred fl ux, which is subject to 

model assumptions ( Comelly_et alj [20 1 ^ _ 

In IFinoguenov et al.l (l2007l 1201 Ol) : Bielbv et aP (1201 Ol) we 
introduced a system of flagging the source identification. The ob¬ 
jects with flag=l are of best quality, with centroids derived from 
the X-ray emission and spectroscopic confirmation of the red- 
shift; flag=2 objects have large uncertainties in the X-ray center 
(low statistics or source confusion) with their centroids and flux 
extraction apertures positioned on the associated galaxy concen¬ 
tration with spectroscopic confirmation; objects with flag=3 still 
require spectroscopic confirmation; objects with flag=4 have 
more than one counterpart along the line of sight; objects with 
flag=5 have doubtful identifications and are only used to access 
systematic errors in the statistical analysis associated with source 

identification. _ 

Using a catalog of iMiller et al.l (1201 3l) we And a number of 
complex radio sources inside the X-ray galaxy groups, the cor¬ 
respondent group ids are: 3, 12, 19 (contains a Wide Angular 
Tail source), 26, 43, 52, 57. All these sources do not have a two- 
dimensional match of the shape of their X-ray emission with the 
radio. In all cases, but group 43, we can also rule out a substantial 
(> 10%) contribution of the IC emission associated with radio 
source to the X-ray flux. For group 43 this contribution can be 
up to 50%, estimated using the part of the source flux in the area 
overlapping with the radio emission. We note that the associated 
with group 43 radio galaxy is the strongest FRII source in CDFS. 
Other studies typically And one IC X-ray source per square de¬ 
gree (iJelic et alJl20I2l) . so the statistics of CDFS is consistent 
with that. 


5.2. Statistical properties of the groups 

In FiglTT]we plot the sample in the flux-redshift plane. The con¬ 
fusion of sources and our approach to reduce it using 32" spatial 
scales for the flux extraction, results in large flux corrections at 
low-z. The correction approaches unity (thus no correction at 
all) for z > 0.5 sources with a high significance of the detec¬ 
tion. This also introduces a redshift dependence to the flux limit, 
with a limiting flux of 10“'^ ergs s“' cm“^ at z=0.05 levelling 
off at 1.5 X 10“'® ergs s“* cm“^ at the redshifts exceeding 0.5. 
However, to account for this effect is straight-forward. Our ex¬ 
perience shows that different science goals require different sub¬ 
samples, a mass-limited sample, for example, would be selected 
differently. Also, some definition of galaxy groups would make 
a cut on X-ray luminosity, removing the need for an equal flux. 
For most of our own work, the high-z galaxy groups are the ones 
that we are most interested in dZinaro et al.ll2013ll2014h . 

The X-ray detected groups span a large range of X-ray lu¬ 
minosities (lO'" - 10"'^ ergs s“'). The total masses of the X- 
ray groups are derived by applying the empiri cal Ly-M?nn re¬ 
lation determined for the COSMOS groups in iLeauthaud et al.l 
(1201 oh via the weak lensing analysis. FiglTS] shows the de¬ 
rived mass range and compares it to t he calibrated range in the 
COSMOS survey (iGeorge et al.ll201 ih . The ECDF-S groups oc¬ 
cupy a unique mass-redshift space, which influence our under¬ 
standing of galaxy evolution in the group environment. This is 
explored in the ded i cated follow-up papers ( Po pesso et al.ll20l3: 
IZinaro et al.l 1201 3l I20l4 lErfanianfar et al.l2014h . The result- 
ing ECDE-S sample of X-ray detected groups ranges between 
5x10'^ and 5 x IO'^Mq. Eor the first time, the derived masses 



Fig. 15: Comparison of the mass-redshift sampling of the ECDF-S 
(filled black circles) and COSMOS (filled grey circles) X-ray 
group samples. Definition of mass with with respect to the criti¬ 
cal density. The ECDF-S groups extend to much lower masses, 
while occupying a similar redshift range. An improvement in 
the mass sensitivity of the survey scales as exposure to the 
power of 3.3, so 30 times deeper data in ECDF-S results in 
a 3 times better mass limit. 


cross the IO'^Mq mass range, much below the typical X-ray 
group mass of 5 x 10*^. 


5.3. Consistency with Cosmology 


In Fig[T3]we compare the masses and redshifts of the detected 
groups with the density of groups, expecte d from the Planck cos¬ 
mology jHanckColISoration|et^|20T3) and the scaling rela¬ 
tions of iLeauthaud et ah ( 20101) . One can see that the two bright 
low-z groups are unusual for the size of the field, while there is 
a lack of structure at 0.2 < z < 0.5. 


Most previous studies, which rep orted the count s from 
extended source s in deep surveys (IG acconi et akl l2002t 
iBauer et aTll2004 IFinoguenov et al.l 1200711201 Ol) primarily re¬ 
port the emission identified with galaxy groups. Also the mod¬ 
elling of logN-logS of extended sources assumes that it stems 
from groups and clusters of galaxies. 

In FiglT7]we show the log(N > S)- log(S ) of X-ray groups 
in ECDF-S. The data are c onsistent with the p rediction of no 
evolution in the XLF from iRosati et al.l (l2002l) down to 10“'® 
ergs s ' cm“^ fluxes, where the predicted number of groups 
is 500 groups per square degrees and the measured values are 
bounded by the 300-700 range. A power law approximation to 
the logN-logS gives an index of -0.85 (or 1.85, conventionally 
used for AGN differential logN-logS). We have not corrected for 
the faint low-z groups that cannot be detected in our survey, but 
this correction is small due to the low volume at low-z. 


We find that the observed counts are consistent with num- 
ber counts predicted for a flat ACDM Planck cosmology 
(IPlanck Collaboration et al.ll2014l) with = 0.3 and h - 0.7 
and cTg = 0.81, when the results of the simulations of the source 
detection in ECDF-S and the scaling relations of Leauthaud et al 
(2010) are combined. We note that the differences in the cosmo¬ 
logical parameters affect onl y mildly the deriva tion of the scal¬ 
ing relations. As explored in iTavlor et al.l (1201 2l) . the sensitivity 
of lensing geometry for COSMOS group experiment to Oa is 
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Fig. 16: The X-ray luminosity function of ECDF-S groups. Black dots 
show the measurement using the full field, while black crosses 
show the measurements excluding the central region where 
there is a low spatial density of groups. Gray crosses show 
the results from the COSMOS field. Dashed and dotted curves 
show the local XLF in the Northern and Southern Hemisphere, 
revealing an effect of sample variance, caused by small vol¬ 
umes probed by RASS at low group luminosities. 


0.15 at 68% confidence level, while the differences to Planck 
cosmology are much smaller, 0.04. 

The predicted number of sou rces in the Planck cosmology 
( Planck Collaboration et al.|l2014h and the scaling relations of 
iLeauthaud et al. ~ if2010h . combined with the presented detailed 
simulations of the source detection in ECDF-S is marginally in¬ 
consistent with the data. Introducing the 30% deviations in the 
scaling relations, allowed by our calibrations, is required to re¬ 
produce the best fit logN-logS. 

FigHSl compares the X-ray luminosity funct ion in ECDF- 
S with that of COSMOS (iFinoguenov et al. 200^_^nd the local 
measurements based on RASS ( Bohringer et al. 2001h . In com¬ 
puting the X-ray luminosity function (XLF), we limit the sample 
to z < 1.2, where our spectroscopic follow-up is complete, and 
we can account for our redshift cut through the volume calcu¬ 
lation. We illustrate the sample variance within the ECDF-S by 
using the full and partial areas of the survey, which also probes 
the importance of the completeness correction. We find the sta¬ 
tistical and systematic errors on the XLF to be similar. We cor¬ 
rect for the detection completeness using the simulations. This 
introduces a different limiting redshift, as a function of luminos¬ 
ity at which the detection is complete. While in the calculation 
of XLF this is simply the effective volume, there is a difference 
in the effective maximum redshift probed by the data as a func¬ 
tion of the luminosity, which limits the statement about the XLF 
redshift dependence. At luminosities near 10"^^ ergs s ', no evo¬ 
lution of XLF betwe en z < 0.6 and 0 . 6 < z < 1.2 has been 
previously shown by IFinoguenov et al.l (l2007h using COSMOS 
data. ECDF-S both extends the measurement of XLF down to 

sampled at z < 0.2 
s ' to a redshift 
of 1.2. So in agreement with the COSMOS data, which sampled 
those systems to a redshift of 0.6, our current work extends the 
claim of no evolution in XLF down to luminosities of 3 x 10^^ 


unprecedented luminosities of 10“^* ergs s * 
and samples groups with of 3 x 10"'^ ergs 



^500 (0-5-2 keV) [ergs cm ® 


Fig. 17: The log(N > S )-log(S ) of X-ray groups. The grey curves show 
ECDF-S data and the Icr envelope shown as dashed curves. The 
solid black curve shows t he prediction of a n on-evolving X-ray 
luminosity function from iRosati et alJ(l2002h . The long-dashed 
line shows the simulated detected cou nts usin g Planck cosmol¬ 
ogy l|Planck Collaboration et al.l2014ll and th dUeauthaud et aP 
( |20W scaling relation. The dotted line, illustrates the effect on 
changing the normalisation of scaling relations by increasing 
the associated mass by 30%, allowed by our calibration at faint 
fluxes (below 10“'^ ergs s“' cm“^). 


ergs s '. We note that this is not a trivial addition to the previ¬ 
ous COSMOS result for 10“^^ ergs s ', given that feedback pro¬ 
cesses are expected to play an important role at low-luminosity 
groups, which might cause differences in the evolution of XLF 
as a function of luminosity. While all dataset probe groups at 
10"'^ ergs s“' luminosity, the maximum redshift for a detecting 
such systems changes from 0.02 for the RASS (and the differ¬ 
ences between North and South can be interpretted as sample 
variance), to 0.3 for COSMOS to 0.6 for the ECDF-S. No evo¬ 
lution at low Lx does not contr adict to the results o n XLF evolu¬ 
tion at Lx > 5 X 10"'"' ergs s ' (iKoens et al.ll2013l) driven by the 
massive cluster growth. 

The conclusion on the absence of strong XLF evolution at 
the luminosities below 10"'^ ergs s“' is in agreement with the 
logN-logS modelling, which is best fit by the non-evolving XLF. 
In the probed range of X-ray luminosities, no detectable evolu¬ 
tion in the XLF is expected from a combination of cosmology 
(reducing the number of groups of a given mass) and evolution of 
scaling relations (increasing the X-ray luminosity of each grou p 
for a given mass) adopted in our work (IFinoguenov et al.ll2010h . 

FigHS] shows the dn/dz/dQ distribution of the ECDF-S 
groups. The grey curve shows the cosmological prediction with 
parame ters fixed to the Planck 1 3 cosmology and the scaling rela¬ 
tions of iLeauthaud et all (1201 Ol) (solid curve) and a 30% change 
in the normalisation of the scaling relations allowed by our cal¬ 
ibration (dashed curve). We conclude that sample variance, dis¬ 
cussed above is caused by the lack of structure at 0.2 < z < 0.5 
and marginally at 1.2 < z < 1.5, while at other redshifts the 
ECDF-S can be considered as a representative field. We further 
note that the modelling of dn/dz/dQ is sensitive to the detection 
of systems at the detection limit. More work on understanding 
the variety of shapes of the intragroup X-ray emission is needed 
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Fig. 18: The dnldzIdO. [deg“^] distribution of the ECDF-S groups 
(black crosses). The prediction fr om the detailed detection 
simul ation and Planck cosmology dPlanck Collaboration et alj 
l2014h is shown as solid grey curve and the effect of 30% change 
in the scaling relations is shown by the dashed grey curve. 


in order to derive conclusions on the cosmological parameters 
implied by the survey. As an example, many of the groups re¬ 
ported here have been previously detected but assigned a much 
smaller flux. On the other hand, some of the new detections have 
fluxes above the formal limits of previous work, illustrating how 
the variety of shape results in the source detectability. While 
we attempt to account for this e ffect, our model parame ters are 
fixed to the local measurements (iFinoguenov et al.ll2007h . which 
might not be representative for the high-z groups. The problem 
with the flux correction is most important for systems at the 
detection limit, as only part of the source is detected. The low 
statistics prevent us from evaluating dn/dz/dQ for the high-flux 
subsample. 


6. Auto-correlation function of groups 

We can use the two-point correlation function to measure the 
spatial clustering of galaxy groups and to estimate their total 
mass. With 40 spectroscopically identified groups we just have 
enough systems to constrain these statistics. We use the same 
random catalog that has been use d throughout the paper and 
apply the Landy-Szalay estimator (iLandv & Szalavl [19931) . To 
separate the effects of redshift distortions we measure the spa¬ 
tial correlation function in projected separations between groups 
in the direction perpendicular (r^) and parallel {n) to the line- 
of-sight. We then integrate over the velocity {n) component of 
the correlation f unction. FigIT9l shows t he projected correlation 
function, Wp(rp) (iDavis & Peeblesiri983l) . which removes the ef¬ 
fect of infall on the clustering signal. In the halo model approach 
the amplitude of the group clustering signal at large scale (two- 
halo term) is related to the typical mass of the galaxy groups 
through the bias factor. In detail, 

Knod(’'p^ = blhs^DM(r^,z=0) ( 1 ) 

In perfor ming this analysis, we merge the groups within r 2 oo 
of each other (lAllevato et al.ll201^ . which removes the one halo 


term in the correlation function. As shown in Fig[T9] at projected 
separations exceeding 1 Mpc h ', the shape of the galaxy groups 
correlation function is well-fit by the two-halo term. The mea¬ 
surement of an excess clustering signal at 0.2 Mpc h“' indicates 
that non-linear gravitational collapse is nevertheless affecting the 
signal, so extension of the comparison between the prediction for 
the linear growth of the two halo term to rp < 1 Mpc h * is not 
supported by the data. 

In modelling we compare the measured amplitude of the 
two halo term with the prediction of linear biasing using 
the mass of each group that contributed to pair statistics and 
weighted with the large-scale structure (LSS) growth function. 
In detail, for each galaxy group at redshift zi, we estimate the 
bias factor corresponding to a DM halo mass M 200 (h * Mpc): 

bi ^ b{M2oo,Zi) ( 2 ) 


where b(M 2 S)o ,Zi) is evaluated fo llowing the bias-mass re¬ 
lation described in ISheth et alJ (1200 ll) . The linear regime of the 
structure formation is verified only at large scales, which is fur¬ 
ther confirmed by our data in Fig[T3 so we estimated the average 
bias of the sample, including only the pairs which contribute to 
the clustering signal at rp - \ - AO Mpc h '. As described in 
lAllevato et alJ (1201 ll) . we define a weighted bias factor of the 
sample as: 


Y,ijbibjDiDj 


'^predicted 


N, 


pair 


(3) 


where bibj (each defined by Eq.2) is the bias factor of the 
pair i-j and Npair is the total number of pairs in the range Vp - 
1-40 Mpc h“'. The D factor is defined by D\{z)ID\{z - 0), 
where D\{z) is the growth function (see eq. 10 in Eisenstein & 
Hu (1999) and references therein) and takes into account that 
the amplitude of the DM two-halo term decreases with increas¬ 
ing redshift. We verified that the bias factor estimated using the 
correlation function of galaxy groups (Eq.l) is consistent with 
the weighted bias factor (bpredicted)- 

Eor all (40) groups with flag < 3 and z < 1.3, we find a 
best fit bias bohs - 2.28 + 0.25 , estimated using a Xc minimi¬ 
sation technique with 1 free parameter, where ^ 

A is a vector composed of Wobsifp) - Mcoi' is the 

covariance matrix. The subscript c denotes that the correlations 
between errors have been taken into account through the inverse 
of the covariance matrix. 

Although the measurement is affected by sample variance, 
we can reproduce its level with the help of Eq.3, as the bias 
prediction is done using the properties of the sample, which 
is at variance with the expectation for an average mass func¬ 
tion and a uniform spatial distribution of groups. We associate 
a mass to each group using the measured X-ray luminos ity of 
groups and the scaling relation of iLeauthaud et al] (l2010l) . The 
error in the prediction is estimated using the scatter in the mass- 
lumi nosity scaling relatio n, constrained from COSMOS to be 
20% (lAllevato et al.ll20l3) . We predict bpredicted - 2.08 + 0.07. 

Although our method can help in the case of small fields, 
any clustering method needs to cover the angular scales corre¬ 
sponding out to projected radii >10 Mpc beyond which the cor¬ 
relation signal drops and the noise estimate is possible. Without 
this the integral constraint affects the measurement. The size of 
ECDE-S is just big enough for such a measurement to succeed. 

In order to verify that our comparison is indeed unaffected 
by LSS, we repeat the analysis excluding groups (3 in total) 
and random objects from the central region, obtaining bobs - 
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Fig. 19: Projected autocorrelation function of X-ray galaxy groups in 
the ECDF-S. The red points show the results for all the groups 
and the black points again for groups but excluding the cen¬ 
tral area from both real and random catalogs. As we discuss 
in the text, we do not see any significant change in the results, 
after introducing the method for correcting the first order ef¬ 
fects from cosmic variance. The dotted lines show the two-halo 
terms b-^redicied^DM(r'p,z = 0) for all groups (black) and exclud¬ 
ing the central area from the real and random group catalogs 
(red). 


2.13 + 0.24 and bpredicted -2.\2± 0.08. Although the agreement 
seems to be better, we point out that within the statistical errors, 
the two measurements are the same. We have also tested our 
method using the Millennium catalogs, comparing the bias of 
subsamples of halos with the predicted bias for the halo masses 
and the bias-mass relation suitable for the cosmological param¬ 
eters of the Millennium simulation, revealing an agreement to 
better than 10%. 

Based on this agreement, we can exclude large (> 30% 
in mass) departures from the scaling relations we use, which 
implies that the ECDF-S sample indeed consists of low-mass 
(IO'^Mq) groups and not of some imaginary low-luminosity 
massive clusters. 


7. Weak Leasing calibration 

In this section, we describe a stacked weak lensing analysis of 
the ECDF-S groups using high-resolution data from the HST 
GEMS survey. 


the PSF and redshift distribution as a function of magnitude and 
briefly summarise here. 

Object cata l ogues are created with SExtractor 
(iBertin & ArnoutsI Il996h and hand masked to remove false 
detections along chip boundaries, star d iffraction spikes, 
satelli te trails, and reflection ghosts as in iMacDonald et al.l 
(l2004l) . Geometric distortions due to the off-axis locatio n of 
ACS are cali brated via a mode l from (jMeureret^ 2()03[) and 
multi drizzle (iKoekemoer et al.l l2003h . Hevmans et al.l (l2005h 
found no evidence of problems arising from charge transfer 
efficienty (CTE), which causes a correlation of object shapes 
with the read-out direction and distance from the read-out 
amplifier, so no correction for CTE is made. However, there 
is a strong anisotropic PSF distortion which must be carefully 
modelled and removed to allow confidence in measured shapes. 

The PSF is characterised through non-saturated point-like 
objects selected via th e stellar locus on the size-magnitude plane. 
iHevmans et al.l (l2005l) fit a two-dimensional second-order poly¬ 
nomial to the anisotropic PSF, modelling each chip and data 
with different depths. The fit is done with a two-step iterative 
procedure with 3cr outlier rejection. After the correction has 
been applied, the residual mean stellar ellipticity is reduced from 
~ 4% to ~ 0.03%, consistent with zero within the error bars. 
G alaxy ellipticities a r e measured using the me t hods described 


in 


Kaiser et al.l (1 19951) i lLuppino & KaisCT(ll997l) : lHoekstra et al.l 
(Il998h (KSB-b) and converted to shear estimates using the pre- 
seeing shear polarizability tensor. The level of shear calibration 
bias from this method has been shown to be ~ 3% on simulations 
(IHevmans et al.ll200^ . which is much smaller than the statisti¬ 
cal uncertainties of this analysis. Source galaxies are selected as 
having size >2.4 pixels, galaxy shear < 1, 24 < meoe < 27, and 
SNR > 15, yielding 41,585 galaxies or a number density of ~ 52 
galaxies per square arcminute. 

Knowledge of the redshift distribution is also important for 
interpretation of the lensing signal. We assume that a magnitude- 
dependent redshift distribution can be parametrised as 


n{z, mag) cx exp 


Zoimag) 


1.5 


(4) 


where zn - Zm/1-412 with Zm being the median redshift 
(iBaugh & Efstathi^ll994l) . ZmO^eoe) is measured for galaxies 
from COMBO-17 with multi-band photometric redshifts and 
galaxies from VVDS with spectroscopic redshifts, and the best 
linear fit is; 

Zm - -3.132 + 0.164m606 


for (21.8 < meoe < 24.4). We extrapolate the above relation¬ 
ship for galaxies fainter than meoe = 24.4, which agrees 
with the Zm-m 606 rel ationship determined for the Hubble Deep 
Field North (HDFN; iLanzetta et al.lll9^lFernandez-Soto^tai, 


[1999). Further details are given in Section 6 of IHevmans et al' 
(I2nn.5h . 


7.1. GEMS Source Catalogues 

The Galaxy Evolution from Morphology and Spect r al energy 
distributions survey (GEMS; iGiavalisco et al.lI2004I : iRix et al.l 
l2004l) consists of deep optical data (5cr point source detection 
limit of OT 606 = 28.3) taken by the Advanced Camera for Surveys 
(ACS) on the Hubble Space Telescope (HST) spanning 795 
arcmin^ c entered on ECDF-S. W e refer the reader to Sections 
3 and 4 of IHevmans et akl (l2005l) for details of the GEMS data 
processing, including the cataloguing, characterisation of both 


7.2. Lensing Signal 
7.2.1. Formalism 

We measure the tangential component of lensing-induced shear 
jT, which is the component of shear perpendicular to the line 
transversely connecting the lens and source positions, yr is re¬ 
lated to the so-called differential surface mass density AS as fol¬ 
lows; 

yr'Scnt = S(< rp) - t(rp) = AS (5) 
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where is the physical transverse separation between the lens 
and source positions, £ is the surface mass density averaged 
within rp, and £(rp) is the mean surface mass density at rp. The 
critical surface mass density Scrit is given by 

c- Ds 
AnG DiDis 

where c is the speed of light, G is the gravitational constant, and 
Di, Ds , and Dis are the angular diameter distances to the lens, 
source, and between lens and source, respectively. 

The weighted lensing signal around each ECDF-S group po¬ 
sition is averaged over bins of rp and can be formally described 
as; 


7t - 


\-ijV|ens \~'^STC 

Hi Hj Wjjjj 
Wjj 


(7) 


Wii = 


(O-SN + 0-?) 


The weights wjj depend on the intrinsic shape noise ctsn and 
the measurement error cTg. The physical scale used to convert 
from angular distances to rp is determined by the spectroscopic 
redshift for the given group lens, given in Tab|4] 

The averaged lensing signal from the ECDF-S groups can 
be co mpared with th e expe cted signal from a dark matter halo 
with a iNavarro et alJ (Il996h density profile. The equations de- 
scribing the radial depen dence of the shear can be found in 


IWright & Brainerd (2000^We fix the concentration using the re¬ 
lation in Duffy et al. (l2008h . effectively turning the NEW model 
into a single parameter profile dependent only on the halo mass. 
In this case, we use M 200 , which is the mass enclosed within a 
sphere with radius R 200 , the radius at which the mean enclosed 
mass density is 200 x pc, and pc is the critical mass density. 


7.2.2. Results 

We measure the mean lensing signal using Equation |7] for the 
ECDF-S groups that have z < 0.8 and flag - 1. The redshift 
limit is chosen because the higher redshift groups have very few 
background source galaxies and thus mostly contribute noise. 
The choice of the flag is to include only those groups with secure 
X-ray centers, as miscent ering issues can ad ditionally bias the 
lensing measurement low (iGeorge et al.l2012l) . We do not further 
address miscentering due to the large statistical errors, while we 
exclude the shear signal below 0.1 Mpc /r '. 

We measure the mean lensing signal using Equation Q for 
the ECDF-S groups that have Zi < 0.7, as the higher red¬ 
shift groups have very few background source galaxies and thus 
mostly contribute noise. Fig. l20l shows yj as a function of dis¬ 
tance to the group center. The errors are given by bootstrapping 
with 1000 resamples. A common systematics test: the 45-degree 
rotated shear is also plotted, and is consistent with zero. We fit 
an NEW profile for a mean zl - 0.7 and Zs = 1.08 using least 
squares optimisation to the tangential shear signal and obtained a 
best fit mass of M 200 = 1.52 x IO^^Mq. The corresponding shear 
profile is overplotted in Fig. |20] The fractional error on the to¬ 
tal mass is 50%. The average mass of the X-ray groups entering 
the stack (IDs: 3, 8, 9, 10, 12, 16, 17, 18, 21, 25, 27, 28, 29, 33, 
34, 35, 39, 44, 48, 49, 50, 61, 68) based on the extrapolation of 
Lx - M 20 Q relation is 1.88 x IO'^Mq. Thus, the mass calibration 
is confirmed through the weak lensing analysis, yet with large 
statistical uncertainty. 



rp (Mpc h 


Fig. 20: Tangential shear signal as a function of radial separation. 
Errors are measured using 1000 bootstrap resamples. 


8. Kurk superstructure 

CDF-S hosts over-densities of gala xies at z = 1.6 id entified in 
the CMOS spectroscopic campaign (iKurk et al.l2009h . Based on 
the deep X-ray data, we find that there are no obvious X-ray 
counterparts of them at a significant S/N. In FiglST] we present a 
detailed map around t he over-den s e regio n. All 5 putative peaks 
in the photo-z map of iKurk et al.l (l2009h are located within the 
area of positive X-ray flux with low S/N. But, as can be seen 
in the map, the X-ray emission from the over-densities, if any, 
would have been confused with nearby X-ray groups. While 
the detection of the sources is confused with foreground struc¬ 
ture, it is still possible to perform the aperture fluxes, placing 
30" aperture on each component. These sources have fluxes 
at or below 2 x 10“*® ergs s“' cm“^. We include the proper¬ 
ties of these sources in the main group catalog. As can be seen 
there, the putative groups typically have M 200 ~ 2 x 10*^4/©. 
We note that the mos t massive X-ray sele cted system in the 
z = 1.6 structure is the lTanaka et alJ (1201 3ah group with a mass 
of M 200 ~ 3 X IO'^Mq. Given the spatial clustering of galaxies, 
the Kurk over-densities are likely a galaxy cluster in formation. 


9. Discussion and conciusion 

We have presented the detection, identification and analysis of 
the extended sources in the deepest X-ray survey to date - the 
ECDF-S. After the careful subtraction of point-like sources in 
the XMM and Chandra data, we extract extended sources in the 
combined X-ray data. The optical counterparts of these sources 
are searched for using the red sequence technique with the deep, 
multi-wavelength data available in the field. A large combined 
effort of spectroscopic follow-up observations allowed us to de¬ 
rive spectroscopic redshifts of a large fraction of the systems. 

The group catalog contains low-luminosity groups that can 
only be found in deep X-ray surveys. By means of stacked weak 
lensing as well as clustering, we have confirmed that these low- 
luminosity systems are indeed low-mass systems. 

According to the hierarchical model of structure forma¬ 
tion, massive galaxies spend most of their lifetime in group¬ 
sized halos (e.g. De Lucia et al. 2012) where environmental 
processes can strongly affect their evolution. Moreover, in the 
local Universe, groups represent the most common environ- 
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Fig. 21: Signal-to-noise map of X-ray residuals in the 0.5-2 keV band, 
after removing the point sources. The black circles show the lo¬ 
cation of the primary identifications. The solid red circle is the 
center of the Kurk superstructure and the solid green sources 
are other peaks on the photo-z at z=1.6 with sufficient spectra 
and the blue circle is a photo-z peak with just one spectrum, all 
found in the map of Kurk et al. Small solid circles indicate the 
spectroscopic members of the 1.6 wall and the dashed red circle 
is the highest significance X-ray source at z=1.6 presented in 
Tanaka et al. (2013). The radii of the circles correspond to 1?200. 
also giving a sense of the angular/physical scale (e.g. Tanaka’s 
group has a radii of 0.7'), with the coordinates of the circle 
centers listed in the Table |4] 


ment of galaxies (Geller & Huchra 1983; Eke et al. 2005). 
Therefore, studying galaxy groups at different cosmic times is 
vital to understand how the environment affects galaxy proper¬ 
ties (e.g. star formation activity and morphology). ECDF-S cat¬ 
alog has already been exploited for the galaxy evolution as wel l 
as AGN stu dies 1^ yerman et al.l (120101) : iPonesso et al.l (l2012h : 
IZinaroet al.l (120131120l4 . 

The CDFS group catalogue presented in this work has al¬ 
ready been crucial to underline the import ance of groups with 
respect to other environments. For example. IZiparo et al.l(l2013h 
show that, opposite to what it is expected for clusters, groups 
lack of any radial trend in galaxy star formation. Nevertheless, 
the star formation activity in galaxy groups is globally sup- 
pressed with respe ct to group-like density regions and the field 
dZiparo et al.l2014l) . suggesting that processes related to a group¬ 
sized dark matter halos are more efficient in qu enching star 
forma tion than purely density related processes. IZiparo et al.l 
(l2014l) also show that X-ray detected groups exhibit the fastest 
evolution in star formation activity, confirming the key role of 
pre-processing (Zabludoff 8c Mulchaey 1998) in the cosmic de¬ 
cline of star formation. 

The properties of the identified groups such as mass-redshift 
distribution are broadly consistent with the Plankl3 cosmology. 
There is a lack of structure in the ECDF-S at 0.2 < z < 0.5, while 
at other redshifts ECDF-S can be considered as a representative 
field. The field can be well described by the non-evolving XLF, 
which predicts 500 groups per square degree at 10“*® ergs s“* 
cm^ flux limit. 

Our successful extended source detections in the ECDF- 
S paves the way for future large area X-ray missions, such as 
Athena and WFXT, providing a realistic input for the modelling 
of source detection. Our experience shows that the point source 


removal is not the major bottleneck and moderate spatial reso¬ 
lution (of e.g. 5") will be sufficient for that. The most impor¬ 
tant issue, however, is precise modelling of the unresolved back¬ 
ground since the flux of extended sources is typically 10% of the 
background. Furthermore, faint sources are often confused by 
the outskirts of nearby extended sources. The removal of these 
extended fluxes is also important from the point of view of th e 
unresolved Cosmic X-ray Background (ICappelluti et al.ll2012l) . 
as it can otherwise be mistaken for the clustering signal of the 
WHIM. The precise modelling of the background will be a ma¬ 
jor challenge in the next generation X-ray surveys. 
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Table 1. Bands employed in the red-sequence technique. 


redshift red-sequence band 

range MUSYC MUSIC 


0.0 < z < 0.2 
0.2 < z < 0.4 
0.4 < z < 0.6 
0.6 < z < 1.0 
1.0 < z < 1.4 
1.4 < z < 2.0 
2.0 < z < 3.0 


U -VysV 
B - Rws R 
V-IysI 
R- z 
I-JvsJ 
z - J ys J 
J - K vs K 


U-V vs V 
B - Rvs R 
B-I vs 1 
R- zvs z 
I-JvsJ 
z - J vs J 
J - K vs K 
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Table 2. Summary of simulations of X-ray group detection in ECDF-S. 


Type 

N images 

N detections 

Sensitivity only 

12972 

9059 

Confusion 

12972 

7561 

Confusion-l-PSF 

12972 

7716 

Scaling 

10262 

7002 

Confusion-l-Scaling 

10262 

5837 


Note. — Each row of the table considers differ¬ 
ent effects. In addition to a full simulation of de¬ 
tection efficiency in the presence of instrumental 
background, foreground and point sources (tagged 
as "Sensitivity only”), we simulate the effect of 
source confusion (tagged as "Confusion”), effect of 
XMM PSF on increased detection due to more flux 
seen in group outskirts ("Confusion-l-PSF”), and 
30% change in the normalisation of Lx - M scal¬ 
ing relation with ("Confusion-rScaling”) and with¬ 
out ("Scaling”) the effect of confusion. 


Table 3. Contamination and completeness of X-ray group detection in ECDF-S. 


flux r < 32" 

contamination 

completeness at z=0.6 

flux significance at z=0.6 

2.e-16 

6% 

50 

3. 

4.e-16 

2% 

90 

6. 

8.e-16 

0% 

100 

12. 


Note. — The important parameter for detection is surface brightness, which is con¬ 
verted to flux using a fixed detection aperture of 32" in radius. The contamination is 
also specific to the detection scales. The simulations of point source confusion are based 
on deeper Chandra data and use the actual XMM ECDF-S mosaic. Completeness is 
estimated based on simulations of halo detection, which exhibits a strong redshift de¬ 
pendence at low redshifts, with milder dependence seen above redshift of 0.6, selected 
for quotation here. Flux significance is quoted for the central part of ECDF-S. 
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Table 4. Catalog of the ECDF-S X-ray selected galaxy groups. See §5.1 for column description. 


O 


ID 

(1) 

lAU Name 

Cl 

(2) 

R.A Decl. 

Eq.2000 
(3) (4) 

z 

(5) 

flux 10“'^ 
ergs s“* 

(6) 

net 

counts 

(7) 

Lo.l-2.4keV 

10'*^ ergs s“‘ 
(8) 

M200 

10*3 Mo 
(9) 

^200 

f 

(10) 

flag 

(11) 

N(z) 

(12) 

^disp 

km s“* 
(13) 

1 

J033225.9-273057 

53.1081 

-27.5159 

0.170 

69.94 ± 4.45 

1666 

8.7 ± 0.6 

7.2±0.3 

4.7 

2 

6 

358 

2 

J033200.4-273444 

53.0017 

-27.5790 

1.470 

4.35 ± 0.63 

404 

97±14 

10±1 

1.1 

3 

3 

515 

3 

J033227.3-274114 

53.1137 

-27.6872 

0.735 

1.36 ±0.22 

2270 

6.3 ± 1.0 

3.5±0.4 

1.2 

1 

17 

313 

4 

J033310.5-274320 

53.2939 

-27.7222 

0.147 

1.24 ±0.27 

478 

0.13 ±0.03 

0.49±0.07 

2.2 

1 

14 

146 

6 

J033156.9-275047 

52.9871 

-27.8463 

0.758 

0.52 ± 0.08 

814 

3.2 ±0.51 

2.2±0.2 

1.0 

5 

3 

270 

7 

J033209.9-274634 

53.0413 

-27.7761 

1.600 

0.31 ±0.06 

587 

18 ±3 

3.2±0.4 

0.7 

1 

3 

354 

8 

J033241.0-274702 

53.1707 

-27.7838 

0.621 

0.12 ±0.07 

218 

0.54 ± 0.32 

0.82±0.28 

0.8 

1 

7 

189 

9 

J033216.0-274944 

53.0669 

-27.8289 

0.667 

0.28 ± 0.07 

579 

1.4 ±0.3 

1.4±0.2 

0.9 

1 

11 

228 

10 

J033200.1-275435 

53.0005 

-27.9097 

0.736 

0.66 ± 0.08 

903 

3.6 ± 0.5 

2.5±0.2 

1.1 

1 

11 

277 

12 

J033256.1-280218 

53.2339 

-28.0382 

0.667 

4.79 ± 0.63 

510 

16±2 

6.7±0.6 

1.6 

1 

10 

384 

15 

J033244.2-273400 

53.1842 

-27.5668 

0.221 

0.64 ± 0.34 

76 

0.17 ±0.09 

0.56±0.17 

1.6 

3 

1 

154 

16 

J033237.3-273526 

53.1556 

-27.5906 

0.188 

2.20 ± 0.32 

456 

0.36 ± 0.05 

0.93±0.08 

2.2 

1 

5 

182 

17 

J033244.1-273928 

53.1839 

-27.6577 

0.149 

1.54 ±0.24 

612 

0.16 ±0.03 

0.57±0.06 

2.3 

1 

25 

154 

18 

J033334.8-273950 

53.3951 

-27.6639 

0.520 

1.56 ±0.48 

114 

2.9 ± 0.9 

2.6±0.5 

1.4 

1 

5 

273 

19 

J033151.4-273952 

52.9643 

-27.6645 

1.036 

0.57 ± 0.10 

447 

8.5 ± 1.5 

3.3±0.4 

0.9 

3 

1 

323 

20 

J033226.4-274031 

53.1101 

-27.6754 

1.041 

1.05 ±0.17 

1070 

13±2 

4.3±0.4 

1.0 

1 

10 

353 

21 

J033246.0-274118 

53.1917 

-27.6882 

0.732 

0.94 ± 0.09 

1318 

4.6 ± 0.5 

2.9±0.2 

1.1 

1 

24 

293 

22 

J033321.4-274124 

53.3393 

-27.6900 

1.151 

0.42 ±0.10 

195 

9.3 ± 2.3 

3.1±0.5 

0.9 

3 

1 

324 

24 

J033157.2-274228 

52.9884 

-27.7079 

0.666 

0.34 ± 0.08 

356 

1.6 ±0.4 

1.5±0.2 

1.0 

2 

9 

235 

25 

J033209.6-274242 

53.0401 

-27.7117 

0.735 

0.61 ± 0.08 

896 

3.2 ± 0.6 

2.3±0.2 

1.0 

1 

15 

272 

26 

J033229.4-274408 

53.1226 

-27.7356 

0.076 

1.87 ±0.53 

319 

0.05± 0.01 

0.28±0.05 

3.3 

4 

3 

120 

27 

J033252.7-274432 

53.2198 

-27.7421 

0.534 

0.37 ± 0.06 

546 

0.86 ±0.15 

1.2±0.1 

1.0 

1 

10 

210 

28 

J033321.6-274836 

53.3399 

-27.8101 

0.127 

27.14 ±0.46 

18135 

1.74 ±0.03 

2.65±0.03 

4.3 

1 

49 

255 

29 

J033218.6-274733 

53.0775 

-27.7924 

0.735 

0.38 ± 0.05 

772 

2.3 ± 0.3 

1.86±0.15 

1.0 

1 

25 

252 

30 

J033150.6-274917 

52.9607 

-27.8215 

0.679 

1.84 ±0.10 

2631 

6.8 ± 0.4 

3.9±0.1 

1.3 

2 

22 

321 

31 

J033151.0-275038 

52.9626 

-27.8440 

0.679 

0.55 ± 0.07 

645 

2.4 ± 0.3 

2.0±0.2 

1.0 

2 

21 

257 

33 

J033223.2-274943 

53.0968 

-27.8285 

0.578 

0.12 ±0.04 

188 

0.46 ±0.15 

0.77±0.15 

0.8 

1 

5 

183 

34 

J033300.9-275023 

53.2536 

-27.8396 

0.128 

2.02 ± 0.25 

788 

0.14 ±0.02 

0.53±0.04 

2.6 

1 

23 

149 

35 

J033313.1-275039 

53.3047 

-27.8441 

0.127 

3.36 ± 0.33 

642 

0.24 ± 0.02 

0.74±0.05 

2.9 

1 

45 

167 

37 

J033316.1-275158 

53.3169 

-27.8661 

0.880 

0.63 ± 0.08 

509 

5.7 ± 0.7 

2.9±0.2 

1.0 

5 

5 

301 

39 

J033305.2-275209 

53.2715 

-27.8692 

0.518 

0.36 ± 0.08 

334 

0.80 ±0.17 

1.15±0.15 

1.0 

1 

1 

207 

41 

J033218.2-275226 

53.0758 

-27.8738 

1.098 

0.10 ±0.04 

131 

3.1 ± 1.2 

1.6±0.4 

0.7 

2 

8 

258 

42 

J033136.7-275233 

52.9028 

-27.8759 

1.050 

0.22 ± 0.08 

143 

4.6 ± 1.6 

2.2±0.5 

0.8 

1 

5 

283 

43 

J033218.6-275415 

53.0775 

-27.9042 

0.965 

0.18 ±0.06 

306 

3.1 ± 1.0 

1.8±0.4 

0.8 

1 

11 

263 

44 

J033205.6-275452 

53.0234 

-27.9146 

0.684 

0.81 ± 0.07 

1178 

3.4 ±0.3 

2.49±0.14 

1.1 

1 

19 

276 

48 

J033230.3-275732 

53.1261 

-27.9588 

0.621 

0.12 ±0.08 

98 

0.5 ± 0.4 

0.8±0.3 

0.8 

1 

8 

188 

49 

J033225.0-275844 

53.1042 

-27.9790 

0.126 

7.30 ± 0.36 

2618 

0.45 ± 0.02 

1.16±0.03 

3.4 

1 

37 

192 

50 

J033210.7-275925 

53.0445 

-27.9902 

0.680 

1.47 ± 0.11 

1185 

5.5 ±0.4 

3.4±0.2 

1.2 

1 

15 

307 

54 

J033254.5-274521 

53.2270 

-27.7557 

1.600 

0.19 ±0.04 

283 

13.2 ±2.9 

2.6±0.4 

0.7 

3 

2 

332 

55 

J033310.7-274620 

53.2945 

-27.7722 

0.522 

1.18 ± 0.09 

1476 

2.25 ± 0.18 

2.24±0.11 

1.3 

2 

9 

258 

56 

J033323.0-274615 

53.3458 

-27.7710 

0.835 

1.74 ±0.13 

1347 

11.0 ±0.8 

4.6±0.22 

1.2 

2 

8 

349 

61 

J033320.0-274332 

53.3333 

-27.7257 

0.521 

0.64 ±0.11 

331 

1.30 ±0.22 

1.66±0.17 

1.2 

1 

14 

230 
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Table 4 (cont’d) 


ID 

(1) 

lAU Name 

Cl 

(2) 

R.A Deck 

Eq.2000 
(3) (4) 

z 

(5) 

flux 10-‘5 
ergs s“' 

(6) 

net 

counts 

(7) 

Lo.l-2.4keV 

lO'*^ ergs s“' 
(8) 

M200 

10'3 Me 
(9) 

^200 

/ 

(10) 

flag 

(11) 

N(z) 

(12) 

^disp 

km s“* 
(13) 

63 

J033220.5-274436 

53.0854 

-27.7433 

0.524 

0.31 ±0.06 

321 

0.70 ±0.13 

1.06±0.13 

1.0 

2 

4 

201 

68 

J033230.4-275309 

53.1267 

-27.8859 

0.645 

0.07 ± 0.04 

122 

0.42 ± 0.22 

0.68±0.22 

0.8 

1 

4 

178 

76 

J033252.6-274228 

53.2192 

-27.7079 

1.028 

0.26 ± 0.05 

286 

4.6 ± 1.0 

2.2±0.3 

0.8 

1 

5 

283 

79 

J033234.6-274933 

53.1443 

-27.8258 

0.542 

0.12 ± 0.04 

203 

0.37 ±0.13 

0.7±0.1 

0.9 

2 

8 

176 

Kurk-1 

J033229.0-274247 

53.1209 

-27.7130 

1.610 

0.21 ± 0.04 

473 

14.5 ± 2.9 

2.8±0.3 

0.7 

2 

5 

339 

Kurk-2 

J033225.3-274513 

53.1055 

-27.7537 

1.610 

0.19 ±0.04 

408 

13.4 ±2.7 

2.6±0.3 

0.7 

2 

2 

333 

Kurk-3 

J033216.1-274630 

53.0673 

-27.7750 

1.610 

0.16 ±0.03 

364 

11.8 ±2.2 

2.4±0.3 

0.7 

2 

1 

324 

Kurk-4 

J033213.6-274353 

53.0569 

-27.7313 

1.610 

0.18 ±0.03 

359 

12.8 ±2.4 

2.6±0.3 

0.7 

2 

2 

329 

Kurk-5 

J033235.8-274246 

53.1492 

-27.7129 

1.610 

0.13 ±0.04 

256 

10.3 ±3.1 

2.2±0.4 

0.6 

2 

3 
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